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B.  Objective 

Studies  of  light-matter  interactions  in  quantum-confined  systems  and  nanostructures 
have  provided  great  insight  into  diverse  and  fundamental  problems  such  as  many- 
body  interactions  and  entanglement.  In  particular,  optical  spectroscopy  has  proved 
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to  be  a  powerful  tool  for  elucidating  electronic  dynamics  owing  to  its  ability  to  access 
information  that  is  difficult  or  impossible  to  obtain  otherwise.  The  knowledge 
obtained  from  optical  spectroscopic  studies  complement  those  obtained  from  other 
techniques  and  provide  a  more  accurate  picture  of  system  dynamics.  For  example, 
electronic  properties  of  materials  are  mainly  determined  by  transitions  within  one 
energy  band,  which  describe  carrier  transport  in  real  space.  On  the  other  hand, 
optical  properties  are  often  related  to  transitions  between  valence  and  conduction 
bands.  However,  a  strict  separation  between  intraband  and  interband  phenomena  is 
impossible  since  both  are  simultaneously  present  and  critical  for  determining  the 
performance  limit  of  future  opto-electronic  devices  needed  in  DoD  operations  and 
future  commercial  products  required  by  the  society  striving  on  modern  technologies. 

In  the  initial  proposal,  we  outlined  two  specific  projects:  (i)  investigation  of  coherent 
electronic  coupling  in  semiconductor  quantum  dot  clusters;  (ii)  investigation  of  spin 
wave  dynamics  in  ferromagnetic  multilayer  structures.  In  the  course  of  the  project, 
we  have  expanded  the  first  direction  to  include  investigation  of  metallic  nanoparticles 
and  clusters  as  well  as  hybrid  nanostructures  consisted  of  both  metallic  and 
semiconductor  components.  In  addition,  we  have  studied  an  emerging  class  of 
atomically  thin  semiconductors  known  as  the  transition  metal  dichalcogenides.  In 
this  report,  we  summarize  the  key  progress  occurred  during  the  past  6  years.  These 
results  have  been  broadly  disseminated  via  conferences  and  scientific  literatures. 
We  choose  a  few  representative  examples  to  discuss  in  this  technical  report. 

C.  Scientific  Accompiishments 

We  will  focus  on  plasmonic  and  spin  wave  projects  in  this  progress  report  on  which 
we  have  made  the  most  significant  advance. 

C.1  Plasmonic  metamolecules  exhibit  tunable  optical  magnetic  resonance 

During  2013,  we  demonstrated  a  unique  metamolecules  that  exhibit  a  strong 
magnetic  dipole  resonance  that  spectrally  overlapped  with  the  electric  dipole.  Such  a 
metamolecule  can  be  used  as  the  basic  building  block  for  photonic  materials  with 
properties  that  do  not  exist  in  nature.  Our  work  was  published  in  Nature 
Nanotechnology.  We  describe  the  work  in  details  below. 

The  lack  of  symmetry  between  electric  and  magnetic  response  in  the  optical 
frequency  range  is  a  fundamental  consequence  of  the  small  value  of  the  fine- 
structure  constant  and  is  directly  related  to  the  weakness  of  magnetic  effects  in 
optical  materials^’ Properly  tailored  plasmonic  nanoclusters  have  been  proposed  to 
induce  artificial  optical  magnetism^'®  based  on  the  principle  that  magnetic  effects  are 
indistinguishable  from  specific  forms  of  spatial  dispersion  of  permittivity  at  optical 
frequencies^.  In  a  different  context,  plasmonic  Fano  resonances  have  raised  great 
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interest,  particularly  for  use  in  sensing  applications  that  benefit  from  sharp  spectral 
features  and  extreme  field  localization®'^^.  In  absence  of  natural  magnetism,  optical 
Fano  resonances  have  so  far  been  based  on  purely  electric  effects.  In  this  work,  we 
demonstrate  that  a  subwavelength  plasmonic  metamolecule  consisting  of  four 
closely  spaced  gold  nanoparticles  (NPs)  support  a  strong  magnetic  response 
coupled  to  a  broad  electric  resonance.  This  finding  represents  an  important  step 
towards  the  quest  for  artificial  magnetism  and  negative  refractive  index 
metamaterials  at  optical  frequencies^®'^®. 
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Figure  1.  Illustration  of  the  scattering  response  of  plasmonic  nanorings,  (a)  Scattering 
cross  section  (SCS)of  a  four-particle  nanoring  system  excited  by  a  time-varying  electric  field 
directed  along  the  side  of  the  square  arrangement  of  particles  (left)  and  by  a  time-varying 
magnetic  field  directed  along  the  normal  to  the  plane  of  the  nanoring  (right).  Colors  indicate  the 
charge  distribution  inside  each  nanoparticle,  (b)  Scattering  response  of  a  symmetric  nanoring  and 
(c)  when  small  asymmetries  are  introduced.  Symmetry-breaking  induces  interference  between 
the  electric  dipole  (p)  and  magnetic  dipole  (m)  modes,  which  gives  rise  to  a  pronounced  Fano 
resonance  in  the  scattering  spectrum. 

The  key  concepts  are  illustrated  in  Fig.  1:  a  perfectly  symmetric  nanoring  (Fig. 
1a)  is  designed  to  support  overlapped  electric  and  magnetic  dipole  resonances  at 
visible  frequencies^^.  When  excited  by  a  local  electric  field,  this  geometry 
supports  a  broad  electric  dipole  response.  When  excited  by  a  local  magnetic 
field,  the  induced  circulating  displacement  current  supports  a  sharper  magnetic 
resonance  centered  in  the  same  frequency  range.  In  an  ideally  symmetric 
configuration,  a  time-varying  electromagnetic  wave  would  excite  both  orthogonal 
modes,  which  additively  contribute  to  the  total  scattering  cross  section  without 
interference  (Fig.  1b).  Flowever,  when  small  asymmetries  are  introduced,  the 
modes  lose  their  orthogonality  and  interfere,  leading  to  a  characteristic  Fano  dip 
in  the  electric  dipole  scattering  signature  due  to  magnetic  resonance,  reflected  in 
the  total  scattering  cross  section  (Fig.  1c). 
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Figure  2.  Experimental  setup  and  scattering  measurements  vs.  simulations  demonstrating 
the  presence  of  a  magnetic-based  Fano  resonance  in  the  asymmetric  nanoring.  Schematic 
representations  of  (a)  AFM  nanomanipulation  and  (b)  optical  scattering  setup,  (c)  and  (d)  AFM 
images  of  the  nanorings  over  1.25|im  by  1.25  pm  area.  Measured  scattering  spectra  (e)  when  the 
four  nanoparticles  are  far  apart  and  (f)  when  arranged  in  a  square  aggregation  with  small  gaps. 
S-polarized  incident  light  beam  has  been  considered  here.  No  analyzer  was  present  along  the 
collection  path,  (g)  and  (h)  calculated  total  scattering  cross  section  corresponding  to  the  geometry 
in  (c)  and  (d),  respectively.  Panel  (h)  compares  calculated  scattering  cross  sections  for  a  perfectly 
symmetric  nanoring  and  an  asymmetric  one. 


We  used  AFM  nanomanipulation  to  assemble  the  nanocluster,  as  illustrated  in 
Fig.  2a.  Initially,  gold  nanoparticles  approximately  100  nm  in  diameter  were 
randomly  dispersed  on  a  substrate.  AFM  was  first  used  to  image  and  identify  the 
positions  of  several  close-by  nanoparticles.  Next,  four  nanoparticles  were  assembled 
into  a  square  geometry.  Scattering  measurements  were  taken  using  a  home-built 
dark  field  microscope  (Fig.  2b).  As  discussed  above,  the  magnetic  mode  can  couple 
to  the  electric  mode  due  to  suitable  asymmetries,  affecting  the  overall  scattering 
signature,  including  the  electric  dipole  component  that  our  measurements  directly 
observe.  The  scattering  spectrum  with  no  polarizer  in  the  collection  path  (Fig.  2e) 
was  dominated  by  a  single,  broad  resonance  centered  at  540  nm,  arising  from  the 
dipole  resonance  of  the  individual  nanoparticles.  When  the  gaps  between  the 
nanoparticles  were  decreased,  a  large  red  shift  of  the  electric  dipole  resonance  was 
observed  (Fig.  2f)  due  to  the  increased  capacitive  coupling  between  neighboring 
nanoparticles.  When  the  gap  size  was  further  reduced  below  10  nm  (AFM  image  in 
Fig.  2d),  a  pronounced  dip  centered  at  740  nm  was  observed  (Fig.  2f).  This  spectral 
feature,  quite  peculiar  for  such  a  subwavelength  structure,  corresponds  to  a  Fano 
resonance,  which  originates  from  the  interference  between  the  radiant  (bright), 
broad  electric  dipole  and  the  sub-radiant  (dark),  narrow  magnetic  dipole  modes  as 
illustrated  in  Fig.  1c.  This  collective  plasmon  response  can  also  be  interpreted  in 
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terms  of  hybridized  states^®,  for  which  the  interference  between  bonding  and 
antibonding  modes  gives  rise  to  the  Fano  resonance. 

In  order  to  confirm  this  intuitive  interpretation  of  the  observed  spectrum,  we 
performed  numerical  simulations  using  the  three-dimensional  finite  element  software 
COMSOL  Multiphysics.  The  computed  total  scattering  cross  sections  are  displayed 
in  Fig.  2g-h  (red  lines),  showing  good  agreement  with  our  measurements  before  and 
after  the  AFM  assembly. 

Our  study  convincingly  demonstrates  that  a  four-particle  nanoring  can  support  a 
strong  optical  magnetic  resonance  overlapped  with  a  broad  electric  resonance.  This 
specific  geometry  avoids  radiative  damping  associated  with  higher-order  scattering 
modes  and  boosts  the  magnetic  response,  despite  the  overall  subwavelength  size  of 
the  nanoring.  In  addition,  the  near-field  interaction  between  the  magnetic  and  electric 
modes,  enabled  by  the  small  asymmetries  in  our  geometry,  has  led  to  the  first 
observation  of  a  magnetic-based  optical  Fano  resonance  in  the  total  scattering 
spectrum  in  a  direction  along  which  the  direct  magnetic  response  would  otherwise 
be  zero.  Metamaterials  with  engineered  electric  and  magnetic  resonances  and  their 
controlled  coupling  may  enable  many  fascinating  applications  in  nanophotonics 
including  cloaking,  high  resolution  imaging,  sensing  and  enhanced  nonlinear  optical 
response. 

C.2  plasmonic  orientation  sensor 

We  have  successfully  demonstrated  a  new  type  of  sensors  that  take  advantage  of 
near-field  coupling  to  activate  a  dark  mode  in  plasmonic  NPs.  This  accomplishment 
resulted  from  the  close  collaboration  between  Prof.  Li  and  Prof.  Shvets’  groups.  We 
reported  this  result  in  a  manuscript  just  published  in  Nature  Photonics  and  two 
conferences. 

A  metallic  nanoparticle  (MNP)  can  enormously  enhance  the  absorption  or  scattering 
cross-section  of  a  barely  visible  emitting/scattering  object  (ESO)  such  as  a  molecule 
or  a  quantum  dot  placed  in  its  close  proximity.  The  near-field  coupling  between  the 
ESO  and  the  MNP  underlying  this  enhancement  has  been  exploited  in  applications 
ranging  from  surface  enhanced  Raman  scattering^^’  and  fluorescence^^’  to 
plasmonic  sensors^^’  and  nanolasers^^'^^.  It  can  also  drastically  modify  the 
extinction  spectrum  of  the  resulting  hybridized  ESO/MNP  nano-system  and  lead  to 
Fano  interference®’ 
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Here,  we  introduce  a  novel  concept,  a  plasmonic  nano-protractor  based  on 
polarization  spectro-tomography  (PST)  that  exploits  Fano  interference  to  achieve 
three-dimensional  super-resolution  of  an  ESO  placed  in  near-contact  with  an 
MNP.  By  utilizing  a  weakly  scattering  plasmonic  nanorod  as  a  “dark”  ESO  placed  in 
close  proximity  with  a  “bright”  MNP,  we  experimentally  measure  its  three- 
dimensional  position  and  orientation  with  respect  to  the  MNP,  as  illustrated  in  Fig.  3. 
The  nano-protractor  can  lead  to  breakthroughs  in  diffraction-unlimited,  in  situ 
imaging  of  complex  molecules  labeled  by  multiple  ESOs.  The  idea  of  transducing 
changes  in  spatial  arrangement  between  MNPs  into  spectral  differences  has  been 
exploited  in  plasmonic  rulers  capable  of  measuring  length  variations  in 
biomolecules^"^’  and  nanoclusters^®.  Here  we  coin  the  term  “plasmonic  nano¬ 
protractor”  to  highlight  the  new  capability  of  our  approach  to  pinpoint  the  spatial 
position  and  orientation  of  the  ESO  relative  to  the  MNP.  Our  technique  may  be 
applied  for  studying  in  situ  conformational  changes  such  as  rotation  and  bending  at 
the  single  molecule  level.  Because  multiple  ESOs  can  have  different  resonant 
frequencies,  this  technique  can  be  readily  extended  to  super-resolution  sensing  of 
large  molecules  labeled  by  several  distinct  ESOs. 


^a)  Minimal  Far-field  Scattering 


Light  4  ESO  Only 


Lighf-F MNPS  ESO 


Figure  3:  Conceptual  schematics  of  the  polarization  spectro-tomography  (PST)  of  a  hybrid 
nano-object  comprised  of  the  “invisible”  emitting/scattering  object  (ESO)  and  a  large  metallic 
nanoparticle  (MNP),  (a)  The  target  ESO  alone  cannot  be  detected  by  light  scattered  in  the  far 
field  due  to  its  small  scattering  cross  section.  When  coupled  to  a  probe  MNP  in  the  near-filed, 
the  MNP/ESO  hybrid  produces  strong  scattered  signal  with  characteristic  asymmetric  Fano 
lineshape.  (b)  Schematic  of  the  proposed  PST:  the  Fano  axis  of  the  MNP/ESO  hybrid  is 
uniquely  determined  by  measuring  the  scattering  spectrum  for  a  range  of  polarizer  (analyzer) 
angles,  9  (</>)  .  (c)  The  strength  of  the  Fano  resonance  is  strongly  dependent  on  the 
polarization  of  the  incident  and  scattered  radiation,  (d)  A  representative  plot  for  a  Fano- 
resonant  polarizability  of  a  MNP/ESO  hybrid  is  used  to  identify  the  orientation  of  the  Fano 
axis  from  the  values  of  (9p,  <t>p). 
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To  mimic  a  barely  visible  ESO,  we  used  a  gold  nanorod  which  possesses  a  narrow 
“dark”  resonance  in  the  visible  wavelength  range  corresponding  to  a  quadrupole-like 
charge  distribution  as  shown  in  Fig.  4a.  While  the  stand-alone  ESO  was  not 
detectable  in  the  spectral  range  of  the  detector  due  to  the  small  scattering  cross 
section  of  the  quadrupole  resonance,  the  ESO/MNP  hybrid  was  easily  detected. 
Experimentally,  the  ESO  (a  gold  nanorod  with  diameter  d  «  22  nm  and  length  L 
«  180  nm)  and  the  MNP  (a  gold  sphere  with  diameter  D«  150  nm)  were  brought 
together  by  pushing  the  nanosphere  toward  the  nanorod  using  atomic  force 
microscope  (AFM)  nanomanipulation^^‘'^°,  as  illustrated  in  Fig.  4b.  Several 
consecutive  steps  of  the  manipulation  process  are  shown  in  Fig.  4c.  As  the  ligand- 
coated  ESO  and  MNP  were  brought  into  near-contact,  the  scattering  spectrum  shown 
in  Fig.  4d  evolved  from  that  of  the  isolated  MNP  (green  line)  to  the  asymmetric  Fano- 
shaped  spectrum  of  the  ESO/MNP  hybrid  (blue  line).  The  Fano  resonance  appeared 
near  ~  700  nm,  which  corresponds  to  the  quadrupole  mode  of  the  nanorod.  The 
precise  alignment  between  the  near-touching  ESO  and  MNP  can  no  longer  be 
resolved  by  AFM  imaging.  Remarkably,  such  information  can  be  obtained  by  applying 
the  PST  technique  to  the  far-field  scattering  experiments  as  we  demonstrate  below. 


Figure  4:  Optical  resonances  and  experimental  assembly  of  a  MNP/ESO  hybrid,  (a)  The  dispersion 
relation  (wavelength  X  versus  propagation  wavenumber  P)  of  the  surface  plasmons  on  an  infinite 
gold  wire  with  a  20  nm  diameter.  Fabry-Perot  resonances  of  the  finite  length  nanorod,  Z,=  180  nm, 
are  defined  by  the  resonance  condition  fi  »  nnjL  where  n=l  (2)  is  the  dipolar  (quadrupolar) 
mode,  (b)  An  illustration  of  the  nanomanipulation  method  using  AFM:  the  gold  nanosphere  (MNP) 
is  brought  in  contact  with  a  gold  nanorod  (ESO).  (c)  AFM  images  of  the  nanosphere  and  nanorod  at 
different  steps  (i-iv)  in  the  assembling  process;  (d)  The  scattering  spectra  of  the  MNP  measured 
before  (green)  and  after  (blue)  the  ESO  was  brought  into  contact.  A  Fano  resonance  emerged  at 
X  «  700  nm  after  the  assembly  was  completed. 
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Figure  5:  Experimental  data  of  PST  on  the  MNP/ESO  hybrid  comprised  of  the  plasmonic  nanosphere  and 
nanorod.  (a)  Schematics  of  the  experimental  set-up  and  the  reference  axes;  (b-d)  Measured  scattering  spectra 
for  varying  values  of  the  polarization  angle,  6,  of  the  incident  light  and  for  varying  values  of  the  analyzer 
angle,  </>,  of  the  scattered  light  with  p-polarized  (0  =  0).  The  sample  was  rotated  180°  around  the  substrate 
normal  between  (b)  and  (d).  (e-g)  Spectral  derivatives  £>1,2(0)  and  £>i(</>)  taken  at  the  wavelength  indicated 
by  the  dashed  line  in  (b-d)  as  a  function  of  the  polarizer/analyzer  angle.  The  shaded  planes  in  the  insets  of  (b- 
d)  are  defined  by  the  Fano  axis  projections  and  the  incident  beams  as  in  (b)  and  (d)  or  by  the  Fano  axis 
projection  and  the  substrate  normal  as  in  (c).  (h)  Determination  of  the  Fano  axis  by  combing  either  Oif  with 
(approach  1)  or  02f  with  </)2f(approach  2  and  substrate). 


The  experimental  set-up  for  the  PST  measurements  and  the  definition  of  the 
polarization  angles  are  illustrated  in  Fig.  5a  and  explained  in  the  caption.  For 
convenience,  0/0  in  the  polarizer/analyzer  reference  planes  were  defined  as  follows: 
0  =  0°  direction  corresponds  to  the  p-polarization  of  the  incident  beam,  and  0  =  0° 
corresponds  to  the  intersection  between  the  incident  plane  and  the  substrate.  The 
total  unpolarized  scattering  spectra /i(l;  0)  =  </(l;  0,0) >0  shown  in  Fig.  5b  were 
measured  as  a  function  of  the  polarizer  angle  0  and  wavelength  1  =  Incjcd.  The 
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second  set  of  spectra  /2(1;  0),  shown  in  Fig.  5d,  correspond  to  the  light  incident  from 
the  opposite  side  of  the  assembled  structures  while  preserving  the  same  incidence 
plane.  This  set  of  data  was  taken  by  rotating  the  sample  by  180°  around  the 
substrate  normal.  From  here  on,  subscripts  1  and  2  refer  to  the  two  sample 
orientations.  Finally,  the  polarized  scattering  spectra,  =  7(1,0  =  0,0)  were 

measured  as  a  function  of  the  analyzer  angle,  0,  using  p-polarized  incident  light  as 
shown  in  Fig.  5c. 

The  magnitude  of  the  Fano  feature  is  quantified  by  the  derivative  of  the  scattering 
spectra  evaluated  at  the  Fano  resonance  wavelength  1^.  The  spectral  derivatives, 
=  ^^i,2(^;^)/^^U=AR3nd  0i(0)  =  0/i(l;0)/01|;i=;i^are  plotted  in  Fig.  5e-5f. 
The  spectral  derivatives  reach  their  maximal  values  when  the  angles  of  the  polarizer 
and  analyzer  are  aligned  with  the  projections  of  the  Fano  axis  in  the  respective 
reference  planes.  By  fitting  harmonic  functions  to  Fig.  5e-5g,  we  extracted  the  angles 
corresponding  to  the  maximal  spectral  derivatives  as  6^p  =  120°  ±  5°,  (pp  =  68°  ±  5°, 
and  02F  =  75°  ±  4°.  Any  two  of  these  three  angles  can  be  used  to  define  the  Fano 
axis. 

As  the  insets  in  Fig.  5b-5d  demonstrates,  the  Fano  axis  direction  can  be 
geometrically  defined  by  the  intersection  of  two  planes.  We  used  two  different 
approaches  to  explicitly  reconstruct  the  Fano  axis.  In  the  first  approach,  we  combined 
01F  and  0F,  where  one  plane  is  defined  by  the  polarization  direction  of  0if  =  120° 
and  the  incident  beam  direction,  and  another  defined  by  the  analyzer  direction  of 
0F  =  68°  and  the  substrate  normal.  The  same  procedure  was  used  to  reconstruct  the 
Fano  axis  direction  using  9^p  and  02f  in  the  second  approach.  We  now  specify  the 
Fano  axis  in  the  lab  frame,  as  defined  in  Fig.  5h,  to  facilitate  direct  comparisons 
between  the  two  approaches.  In  this  lab  frame,  the  Fano  axis  direction  extracted  from 
the  first  (second)  approach  corresponds  to  0  =  111°,  0=  68°  (0  =  115°  and 
O  =  74°)  as  indicated  by  the  cyan  (magenta)  lines  in  Fig.  5h.  Flere  (0,0)  are  the 
standard  spherical  coordinates.  Note  that  the  Fano  axis  extracted  from  these  two 
independent  approaches  agrees  within  the  experimental  errors,  thereby  providing  the 
key  experimental  confirmation  of  the  PST  technique  introduced  in  this  work.  Below 
we  demonstrate  how  the  Fano  axis  direction  can  be  used  to  determine  physically 
relevant  quantities  such  as  the  nanorod/nanosphere  contact  point  and  the  relative 
nanorod  orientation.  In  another  recent  experiment,  we  have  extended  this  plasmonic 
protractor  concept  to  a  system  consisted  of  a  plasmonic  nanosphere  and  a  single 
semiconductor  dot.  In  that  case,  the  Fano  resonance  is  mediated  by  single  photon 
absorption  events.  The  Fano  resonance  is  present  even  when  the  QD  has  already 
photobleached  because  of  its  still  intact  ability  to  absorb  photons. 
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C.3  Progress  on  Plasmonic  circuit  project 

We  experimentally  demonstrated  that  circuit  concepts  can  be  used  to  model 
photonic  components,  extending  the  powerful  concept  of  modular  and  lumped  elements 
from  electronics  to  photonics.  This  work  was  published  in  Nature  Communication  in 
2014. 


I 

<0 

b 

X. 

Figure  6:  Assembly  and  characterization  of  a  modular  optical  nanocircuit. (a)  AFM  and  (b) 

dark-field  scattering  images  of  plasmonic  and  dielectric  NPs  randomly  distributed  on  a  glass 
substrate.  Only  Au  NPs  yield  strong  scattering  signals.  Scale  bar,  3  pm.  The  NPs  enclosed  in 
circles  in  a  are  missing  in  b,  thereby  identifying  them  as  dielectric  NPs.  (c)  Illustration  of  AFM 
nanomanipulation  as  a  way  to  dynamically  assemble  complex  optical  nanocircuits  from 
independent  optical  inductors  (plasmonic  NPs,  yellow)  and  capacitors  (dielectric  NPs,  green),  (d) 
Schematic  representation  of  our  dark-field  scattering  measurements.  Light  impinges  at  an 
incidence  angle  of  60  degrees  and  the  scattering  signal  is  collected  along  the  substrate  normal. 


Here  we  employ  small  colloidal  metallic  and  dielectric  nanoparticles  (NPs),  ideal 
platforms  to  translate  these  concepts  to  nanophotonic  systems.  With  optical  nanocircuit 
concepts"^^'"^"^,  we  are  able  to  control  the  intrinsic  optical  impedance  of  NPs  by  their 
geometry  and  material  composition.  This  impedance  is  defined  as  the  ratio  of  the  local 
potential  difference  V=|E|-L  and  the  flux  of  displacement  current  ld=-/w|E|S  through  the 
NP"^^'"^"^,  where  E  is  the  local  electric  field  vector,  L  is  the  NP  length  along  the  electric 
field,  (jo  is  the  frequency  of  operation,  £  is  the  NP  dielectric  constant  and  S  is  its 
transverse  cross-section.  It  follows  that  a  dielectric  NP  with  Re(£)  >0  behaves  as  a 
nanocapacitor,  while  a  metallic  NP  with  Re(£)  <0  acts  as  a  nanoinductor.  Ohmic  loss  in 
the  materials  takes  the  role  of  a  nanoresistor. 

In  this  experiment,  gold  NPs  (~60nm  in  diameter)  and  AI2O3  NPs  (~45nm  in 
diameter)  were  dispersed  on  a  glass  substrate.  The  AFM  image  of  the  sample  (Fig.  6a) 
allowed  us  to  locate  each  NP  but  did  not  distinguish  between  different  NP  types  without 
ambiguity.  Then  we  mapped  the  same  area  using  optical  dark-field  microscopy,  from 
which  only  scattering  signals  from  Au  NPs  were  detectable  (Fig.  6b).  Correlating  the 
AFM  and  optical  scattering  images  allowed  us  to  differentiate  metallic  NPs 
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(nanoinductors)  from  dielectric  NPs  (nanocapacitors)  with  certainty.  By  measuring  the 
scattering  from  each  NP,  we  detected  differences  in  the  scattering  spectra  because  of 
slight  variations  in  NP  shape  and  size.  We  then  picked  a  few  chosen  NPs  with  the 
desired  impedance  and  assembled  them  to  produce  more  complex  connections  (Fig. 
6c)  using  the  AFM  tip.  Finally,  dark-field  scattering  measurements,  shown  schematically 
in  Fig.  6d,  were  used  to  confirm  the  functionality  of  the  nanoclusters  by  comparing  the 
measured  scattering  spectra  with  the  predicted  nanocircuit  response. 


Figure  7:  Second  and  third  order  iumped  nanofiiters.  (a)  AFM  image  and  geometry  of  a  3-NP 
cluster  composed  of  an  AI2O3  NP  sandwiched  between  two  AuNPs.  All  dimensions  are  in 
nanometres.  Scale  bar,  100  nm.  (b)  Dark-field  scattering  measurements  for  s-polarized  light 
exciting  the  nanocircuit  with  the  electric  field  along  x  (red  line)  or  y  (blue),  (c)  Corresponding  full- 
wave  simulations  with  electric  field  along  x  (red  line)  or  y  (blue).  (d,e)  Simulated  electric 
displacement  field  distributions  and  field  vectors  at  the  resonance  wavelengths  for  the  two 
incidence  directions:  K-532  nm  (d)  and  A=524  nm  (e).  (f,g)  Theveninnanocircuit  models  of  second 
orde  and  third  order  lumped  nanofilters,  (h)  Corresponding  circuit  theory  predictions:  X-circuit  (red 
line)  and  Y-circuit  (blue). 

We  first  consider  a  three-particle  nanocircuit  consisting  of  a  dielectric  AI2O3  NP 
sandwiched  between  two  Au  NPs,  as  shown  in  Fig.  7a.  The  measured  scattering 
spectra  from  this  nanocluster  using  s-polarized  light  are  shown  in  Fig.  7b,  with  the 
incident  electric  field  either  perpendicular  to  the  NP  array  (Y-circuit,  blue  line)  or  along  it 
(X-circuit,  red  line).  Owing  to  the  anisotropy  of  this  nanocluster,  the  optical  response 
exhibits  a  spectral  shift  dependent  on  the  direction  of  the  impinging  field  relative  to  the 
NP  array.  A  red  shift  in  resonant  wavelength  is  expected  for  the  X-circuit  (Fig.  7b), 
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which  agrees  with  full-wave  simulations  (Fig.  7c).  Figure  9d,e  shows  the  corresponding 
simulated  displacement  vector  distributions  at  the  two  resonant  wavelengths:  A=532  nm 
for  incident  electric  field  parallel  to  the  x  axis  (Fig.  7d)  and  A=524nm  for  electric  field 
along  the  y  axis  (Fig.  7e).  Similar  to  a  conventional  electronic  circuit,  the  intrinsic 
nanoinductance  of  the  gold  NPs  is  not  affected  by  other  components  in  the  circuit. 

On  the  basis  of  the  different  electric  displacement  field  distributions  (Fig.  7d,e),  two 
circuits  with  different  connections  among  the  elements  (Fig.  7f,g)  account  for  the 
change  in  the  scattering  spectra.  The  three  nanocircuit  elements  are  connected  in 
series  (X-circuit  in  Fig.  7f)  for  incident  electric  field  polarized  along  the  axis  of  the  NP 
cluster,  since  this  excitation  ensures  the  continuity  of  the  displacement  current  across 
the  cluster,  as  confirmed  by  our  simulations  in  Fig.  7d.  In  conventional  circuit  theory,  the 
number  of  independent  reactive  elements  determines  the  order  of  a  filter.  Accordingly, 
the  X-circuit  (Fig.  7f)  realizes  a  second  order  filter,  which  is  an  LC  circuit  formed  by  the 
series  combination  of  an  inductor  (the  Au  NPs)  and  a  capacitor  (determined  by 
the  AI2O3  NP  and  the  fringing  fields).  Conversely,  for  incident  electric  field  perpendicular 
to  the  axis  of  the  NP  cluster  (Fig.  7e),  each  element  experiences  the  same  potential 
difference  (that  is,  voltage),  leading  to  a  parallel  connection  between  them.  This  Y- 
circuit  (Fig.  7g)  forms  a  third  order  nanofilter,  in  which  the  fringe  capacitance  is 
connected  in  series  to  the  parallel  of  the  Au  NP  inductors  and  the  AI2O3  NP  capacitor. 
The  scattering  spectra  predicted  by  our  circuit  model  (Fig.  7h)  are  indeed  quantitatively 
consistent  with  the  measured  and  simulated  spectra  for  both  X-  and  Y-circuits  (Fig. 
7b, c).  By  realizing  different  circuit  configurations  controlled  by  the  direction  of  the 
excitation  field,  we  are  able  to  experimentally  prove  the  design  and  operation  of  optical 
stereo-circuits,  which  do  not  have  a  counterpart  in  the  electronic  realm. 
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Figure  8:  A  fourth  order  lumped  nanofilter.(a)  AFM  image  and  geometry  of  a  4-NP 
complex  nanocircuit  composed  of  three  Au  NPs  and  one  AI203NP.  All  dimensions  are 
in  nanometres.  Scale  bar,  100  nm.  (b)  Dark-field  scattering  measurements  for 
polarization  rotating  from  s  to  p;  (c)  corresponding  full-wave  simulations.  Curves 
in  b,c  have  been  shifted  vertically  to  facilitate  comparisons.  (d,e)  Simulated  electric 
displacement  field  distributions  and  field  vectors  for  polarizations  parallel  (d)  and 
orthogonal  (e)  to  the  array  axis,  at  the  A=688  nm.  (f,g)  Theveninnanocircuit  models,  (h) 
Corresponding  circuit  theory  predictions  for  the  two  orthogonal  polarizations. 

Then  we  further  add  another  metallic  NP  below  the  previous  nanocircuit  design,  as 
illustrated  in  Fig.  8a.  The  scattering  measurements  shown  in  Fig.  8b  were  taken  with 
incident  light  direction  fixed  while  gradually  changing  the  polarization  from  p  to  s.  Our 
experiments  (Fig.  8b)  agree  well  with  full-wave  simulations  (Fig.  8c).  Shown  in  Fig. 
8d,e  are  the  simulated  displacement  vector  distributions  at  the  longer  wavelength 
resonance  (A=688  nm),  for  polarizations  parallel  (Fig.  8d)  and  orthogonal  (Fig.  8e)  to  the 
array,  highlighting  the  polarization-dependent  interaction  among  NPs.  Even  in  this  more 
complex  configuration,  the  electric  displacement  vector  distributions  fully  support  the 
circuit  models  reported  in  Fig.  8f,g,  which  correspond  to  fourth  order  nano-filters.  This 
example  shows  that  we  can  quantitatively  control  the  spectral  response  by  connecting 
an  additional  nanoinductor.  More  specifically,  we  have  been  able  to  add  a  new  zero 
and/or  pole  to  the  scattering  response  by  adding  an  independent  reactive  element  to  the 
circuit  of  Fig.  7.  Rather  than  modeling  a  distributed  photonic  system  with  an  equivalent 
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circuit'^^''^^,  our  work  shows  that  it  is  actually  possible  to  synthesize  photonic  circuit 
functionality  by  assembling  modular  lumped  elements. 

We  would  like  to  emphasize  the  key  message  of  this  work:  the  nanocircuit 
paradigm  provides  an  unmatched  degree  of  simplicity  to  control  and  tailor  the  optical 
response  of  rather  complex  nanophotonic  structures.  Remarkably,  the  impedance  of 
each  nanosphere  remains  independent  of  the  cluster  geometry  and  surrounding 
environment  in  all  the  circuits  we  have  demonstrated,  starting  from  the  isolated  NP  to 
the  complex  fourth  order  nanocircuit.  In  other  words,  each  NP  carries  inherent  optical 
impedance,  which  is  a  property  of  the  particle  itself,  allowing  us  to  modularize  its 
response  and  combine  it  in  complex  configurations. 


C.4  Progress  on  plasmonic  nanoparticle-quantum  dot  hybrid  structure 

We  experimentally  demonstrate  for  the  first  time  that  a  single  semiconductor 
quantum  dot  placed  in  nanometer-scale  proximity  of  a  plasmonic  cavity  can  be  used  to 
control  the  scattering  spectrum  and  anisotropy  of  the  latter.  This  work  was  published  in 
Proceedings  of  the  National  Academy  of  Sciences  in  2015. 

Many  quantum  network  and  information  processing  schemes  require  the  enhanced 
light-matter  interaction  between  a  single  quantum  emitter  and  a  cavity,  enabling  the 
effective  conversion  between  photonic  and  matter-based  quantum  states'^^'^^  Those 
cavity-quantum  electrodynamics  (QED)  effects  require  a  high  Purcell  factor  Fp  <xQ/v, 
where  Q\s  the  quality  factor,  and  V\s  the  volume  of  the  cavity  mode  The  high  Q  of 
conventional  photonic  cavities  is  required  to  compensate  for  relatively  large  (diffraction- 
limited)  mode  volumes  and  comes  at  a  cost:  the  narrow  linewidth  of  cavity  modes 
places  stringent  requirements  on  their  spectral  alignment  with  the  frequencies  of 
quantum  transitions.  Plasmonic  cavities,  on  the  other  hand,  achieve  high  values  of 
Fpwhile  maintaining  moderate  Q  values  because  of  their  ultra-small  modal  volume.  The 
relaxed  spectral  alignment  requirements  facilitate  the  experimental  realization  of  various 
quantum  phenomena,  such  as  collective  photon  emission  from  a  small  ensemble  of 
emitters  and  single  photon  sources  with  tunable  statistical  properties 

Prior  experiments  exploring  cavity  QED  effects  associated  with  single  emitters 
coupled  to  plasmonic  cavities  or  waveguides  focused  almost  exclusively  on  the 
observations  of  reducing  the  emitter’s  lifetimes^®’  The  realms  of  quantum 

information  science  and  plasmonics  have  also  been  bridged  by  demonstrating  that 
photon  emission  statistics,  such  as  anti-bunching  behavior  in  the  second  order 
correlation  function  for  single  photon  sources,  remain  intact  following  the  photon- 
plasmon-photon  conversion  process  The  possibility  of  controlling  the  scattering  of 
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a  plasmonic  nanocavity  by  a  single  (and  inherently  quantum  and  nonlinear)  two-level 
system  has  also  been  proposed^'^’ but  never  experimentally  observed. 

The  strongly  coupled  MNP-QD  hybrid  structure  is  assembled  into  a  well-controlled 
geometry  using  the  technique  of  AFM  nanomanipulation.  The  strong  coupling  between 
the  MNP  and  QD  is  experimentally  confirmed  by  measuring  the  exciton  lifetime. 
Analyzing  the  polarization  and  spectral  properties  of  light  scattered  by  the  MNP-QD 
hybrid,  we  observe  that  the  overall  plasmonic  cavity  scattering  is  significantly  modified 
over  a  broad  spectral  range.  A  Fano  resonance  spectrally  aligned  with  the  QD’s 
quantized  exciton  resonance  is  clearly  identified  when  the  polarization  of  the  scattered 
photon  is  along  the  Fano  axis^^  connecting  the  MNP’s  center  with  the  QD.  The 
anisotropic  scattering  spectrum  observed  in  our  experiments  suggests  that  a 
polarization-controlled,  versatile  quantum  light  source  may  be  realized  in  this  simple 
QD-MNP  cavity  system. 

The  calculated  polarization-resolved  scattering  spectra  by  the  QD-MNP  (diameters: 
2rQD  =  6nm  and  =  30nm)  hybrid  are  shown  in  Fig. 9(a)  for  three  polarization 
angles  0^  of  the  analyzer  placed  in  the  collection  path  of  the  scattering  signal  to  mimic 
the  experimental  setup.  In  the  absence  of  the  QD,  all  scattering  spectra  from  a  single 
MNP  are  independent  of  0^  and  possess  a  single  broad  peak  atl^^p«  520nm 
corresponding  to  the  plasmonic  dipole  resonance  of  the  MNP.  The  introduction  of  a  QD 
under  the  MNP,  with  the  separation  gap  oi  g  =  Inm,  modifies  the  scattering  spectrum:  a 
sharp  Fano  feature  emerges  at  the  exciton  transition  wavelength  =  550nm.  The 
magnitude  of  the  feature  is  a  strong  function  of  the  analyzer  orientation.  If  the  projection 
of  the  Fano  axis  onto  the  analyzer  plane  is  perpendicular  to  the  analyzer  direction 
(0^=7r/2  in  Fig. 9(a)),  then  no  Fano  feature  is  predicted  by  our  calculation.  The 
strongest  Fano  feature  is  observed  for0^  =  0,  and  a  weaker  but  finite  Fano  feature  is 
observed  for  intermediate  angles.  Therefore,  what  is  originally  an  isotropic  scatterer  (a 
spherical  MNP)  is  transformed  into  a  highly  anisotropic  one  by  the  strong  hybridization 
between  the  QD  and  the  MNP. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Fig  9:  Calculation  demonstrating  how  the  near-field  coupling  modifies  the 
far-field  scattering  spectra  of  a  NP-QD  hybrid  structure  piaced  on  a  glass 
substrate. (a)  The  scattering  spectra  of  a  NP-QD  hybrid  excited  by  the  un¬ 
polarized  evanescent  wave  coming  from  the  glass  substrate  side  in  all 
azimuthal  angles.  The  angle  (p^  indicates  the  orientation  of  the  analyzer  in  the 
path  of  the  scattered  light.  Fano  feature  is  the  most  (least)  prominent  when  the 
orientation  of  the  analyzer  is  parallel  (perpendicular)  with  the  in-plane 
component  of  the  Fano  axis,  which  connects  the  QD  and  MNP  centers,  (b) 
and  (c)  shows  the  field  near  the  NP-QD  hybrid  at  500  nm  and  552  nm, 
respectively,  as  indicated  by  the  black  arrows  in  the  scattering  spectrum  at 
(p^  -  0  \n  (a).  The  scattering  signal  at  these  two  wavelengths  is  the  same 
(indicated  by  the  dotted  line  on  the  blue  curve)  while  the  MNP  is  excited  much 
more  strongly  in  Fig  1c.  This  wavelength  dependence  proves  that  the 
presence  of  the  QD  indeed  controls  the  MNP’s  scattering  and  anisotropy 
resonantly. 


The  ability  of  a  single  QD  to  modify  the  scattering  spectrum  of  a  much  larger  MNP 

possessing  a  scattering  cross  section  of  (riyiNpAqo)  ~  15,000  -fold  greater  in 
magnitude  than  that  of  the  QD  seems  surprising.  Naively,  one  may  expect  the  effect  of 
a  QD  on  the  scattering  spectrum  of  the  hybrid  system  to  be  very  small.  As  pointed  out 
by  a  number  of  theoretical  studies^'^’  and  our  own  numerical  calculation  shown  in 
Fig.  9,  this  expectation  is  not  correct.  While  theexciton  dipole  moment  is  too  small  to 
produce  significant  far-field  scattering  on  its  own,  it  is  sufficient  to  de-polarize  the  MNP 
in  the  near-field,  thereby  dramatically  modifying  the  electric  polarizability  of  the 
combined  hybrid  system  at  the  frequency  of  the  exciton  transition. 
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To  illustrate  this  point,  the  near-field  distributions  were  calculated  for  =  500  nm 
(Fig.  9b)  and  X2  -  552  nm  (Fig.  9c),  respectively.  These  two  wavelengths  were  chosen 
because  the  scattering  intensities  are  the  same.  The  much  higher  (by  almost  a  factor 
of  2)  electric  field  induced  on  the  MNP’s  surface  atl  =  l2>  is  offset  by  strong  near-field 
depolarization  (light-color  area  near  the  QD)  of  the  MNP  by  the  exciton’s  dipole. 
Because  the  electric  field  of  a  dipole  rapidly  decays  with  distance,  such  extreme  de¬ 
polarization  (which  can  be  alternatively  interpreted  as  the  excitation  of  high-order 
multipoles  of  the  MNP  by  an  exciton)  can  only  occur  if  the  QD  is  placed  within 
nanometers  from  the  MNP.  Therefore,  it  is  extremely  crucial  to  precisely  position  the  QD 
near  the  MNP  as  accomplished  in  our  experiments. 

We  assemble  the  hybrid  structure  using  the  technique  of  AFM  nanomanipulation 
The  assembly  process  begins  by  dispersing  MNPs  and  QDs  on  a  glass  substrate 
randomly.  We  then  simultaneously  obtain  an  AFM  topography  image  and  a 
photoluminescence  (PL)  image  by  scanning  the  sample  on  a  home-built  integrated 
AFM-confocal  microscope.  We  locate  isolated  MNPs  and  QDs  in  close  proximity  via  the 
AFM  topography  image.  We  then  manipulate  a  near-by  MNP  to  approach  the  chosen 
QD.  We  then  measure  the  lifetime  to  confirm  that  the  MNP  is  indeed  in  the  close 
proximity  of  the  QD.  Lifetime  is,  in  fact,  a  rather  accurate  way  to  characterize  the 
distance  between  the  MNP  and  QD  as  demonstrated  in  our  previous  work 

Measured  ensemble  absorption  and  PL  emission  spectra  (taken  in  solution)  of  QDs 
are  shown  in  Fig.  10b.  The  absorption  spectrum  features  multiple  discrete  exciton 
resonances  at  lower  energies  and  a  continuous  absorption  spectrum  at  energies  above 
the  band  gap  of  the  crystal.  While  all  absorption  resonances  may  influence  the 
scattering  spectrum  of  the  hybrid  structure,  we  focus  on  the  lowest-energy  exciton  state 
centered  near  «  615  nm  with  an  ensemble-averaged  spectral  full-width  at  half¬ 
maximum  (FWFIM)  nm. 

The  dark-field  scattering  experiments  are  performed  using  a  home-built  optical 
system  optimized  for  small  MNP  measurements  (Fig.  10a).  An  un-polarized  white  light 
source  incident  in  a  conical  geometry  generates  evanescent  fields  and  excites  the 
hybrid  structure  from  all  directions.  An  analyzer  is  placed  in  the  scattered  light’s  path  to 
select  the  polarization  of  the  collected  scattering.  A  series  of  such  spectra  are  displayed 
in  Fig.  10c  as  a  function  of  the  analyzer  angle.  A  very  sharp  Fano  feature  can  be  clearly 
observed  for  0^  =  0i  =  30°.  Weaker  Fano  features  are  observed  for  other  polarizer 
orientations,  all  in  qualitative  agreement  with  our  theoretical  predictions  in  Fig.  9a. 
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Fig.  10:  The  scattering  and  anisotropy  of  a  piasmonic  cavity  (i.e.  a  MNP) 
controiied  by  a  singie  QD.  (a)  Experimental  schematics  for  scattering  measurements; 
(b)  Experimentally  measured  absorption  (black)  and  PL  emission  (red)  spectra  from  an 
ensemble  of  colloidal  CdSe/ZnS  core-shell  QDs.  (c)  Measured  scattering  spectra  of  the 
assembled  QD-MNP  hybrid  structure  at  different  analyzer  angles.  The  Fano  resonance 
indicated  by  the  dotted  vertical  line  spectrally  aligns  with  the  lowest  exciton  state 
measured  in  the  absorption  from  an  ensemble  of  QDs.  The  Fano  resonance  is  most 
pronounced  at  0^  =  30°.  (d)  Measured  scattering  spectra  of  a  bare  30  nm  diameter 
MNP  at  different  analyzer  angles.  The  small  changes  near  550  nm  or  below  are  likely 
due  to  deviations  from  a  perfectly  spherical  shape  of  the  MNP. 


To  further  confirm  that  the  polarization  dependence  in  the  scattering  spectra  of  the 
hybrid  structure  indeed  originates  from  the  coupling  between  the  QD  and  MNP,  we 
show  the  scattering  spectra  of  a  bare  MNP  in  Fig.  lOd.  The  spectra  do  not  display  any 
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Fano  features  or  analyzer  angle  dependence  in  the  spectral  proximity  of  the  exciton 
resonance  Iq/j.  Therefore,  it  is  indeed  the  coupling  between  a  single  QD  and  the  MNP 
that  turns  an  otherwise  isotropic  plasmonic  cavity  into  a  strongly  anisotropic  one. 
Because  a  single  quantum  absorber  achieves  this  effect,  one  can  envision  the  proposed 
hybrid  system  as  an  experimental  platform  for  observing  a  plasmonic  cavity  anisotropy 
controlled  by  optical  nonlinearity  at  the  single-photon  level.  The  small  angular  variations 
of 5(1,0^)  at  shorter  wavelengths  (around  550nm  or  below)  in  Fig.  10c  and  Fig.  lOd 
most  likely  arise  from  a  small  intrinsic  deviation  of  the  MNP’s  shape  from  an  ideal 
sphere.  Unlike  the  extrinsic  anisotropy  induced  by  the  QD,  it  cannot  be  optically 
controlled  and  is  not  of  interest  for  nonlinear  quantum  optics. 

In  summary,  we  have  demonstrated  that  a  single  semiconductor  QD  coupled  to  a 
MNP  cavity  can  effectively  control  the  scattering  spectrum  of  the  latter,  as  well  as  render 
it  highly  anisotropic.  The  speculative  implications  of  such  extrinsic  anisotropy  are  very 
intriguing.  Qn  the  one  hand,  it  serves  as  an  orientation  sensor  to  determine  the  relative 
locations  of  the  QD  and  MNP.  Qn  the  other  hand,  it  should  be  possible  to  observe 
polarization-dependent  photon  statistics  of  light  scattered  from  the  QD-MNP  nano¬ 
hybrid. 


C.5  Re-Examine  Fundamental  Optical  Properties  of  Silver 

In  this  work,  we  explore  the  intrinsic  fundamental  optical  properties  of  Ag.  Using 
atomically  smooth  epitaxial  Ag  films,  we  extracted  new  optical  permittivity  highlighting 
significant  loss  reduction  in  the  visible  frequency  range.  This  work  was  published  in 
Advanced  Materials  in  2014. 

Loss  represents  the  most  serious  challenge  that  impedes  progress  and  broad  impact 
towards  practical  technology  in  the  field  of  plasmonics®®'®®.  Silver  {Ag)  is  by  far  the 
preferred  plasmonic  material  at  optical  frequencies  due  to  its  lowest  loss  among  all 
metals.  However,  large  discrepancies  exist  among  widely  quoted  values  of  optical 
permittivity  in  Ag  due  to  variations  in  sample  preparation.  Here  we  push  to  the  limit  of 
intrinsic  fundamental  optical  properties  of  Ag.  We  extracted  new  optical  permittivity 
highlighting  significant  loss  reduction  in  the  visible  frequency  range  from  atomic  smooth 
silver  film.  We  measured  a  largely  enhanced  propagation  distance  of  surface  plasmon 
polaritons  (SPPs),  which  confirmed  that  the  intrinsic  loss  in  Ag  is  lower  than  previously 
considered  possible.  The  new  optical  constants  are  free  of  extrinsic  spectral  features 
typically  associated  with  grain  boundaries  and  localized  plasmons  inevitably  present  in 
thermally  deposited  films. 

In  this  work,  we  performed  careful  spectroscopic  ellipsometry  (SE)  measurements  and 
analyses  on  atomically  smooth,  epitaxially  grown,  single  crystalline  Ag  films^^'^\  and 
accurately  extracted  Kramers-Kronig  (K-K)  consistent  optical  constants.  Qur 
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measurements  suggest  that  the  intrinsic  loss  in  Ag  is  significantly  lower,  by  up  to  a 
factor  of  2  in  the  visible  wavelength  range,  than  the  best  values  previously  reported  by 
JC^^.  We  also  measured  SPP  propagation  distances  along  these  epitaxial  Ag  films, 
finding  greatly  enhanced  propagation  lengths,  approaching  the  fundamental  limit 
determined  by  the  new  optical  constants  at  both  visible  and  near-infrared  (NIR) 
frequencies.  The  measured  propagation  distance  confirms  that  the  loss  determined 
from  the  permittivity  data  reported  by  JC  is  not  an  intrinsic  limit. 

Atomically  smooth  Ag  films  were  grown  using  molecular  beam  epitaxy  (MBE)  on 
heavily  doped  Si(111)-7x7  substrates.  The  high  surface  quality  of  a  45  nm  oxide  capped 
expitaxial  film  was  confirmed  by  atomic  force  microscopy  (AFM)  shown  in  Fig.  11a.  For 
comparison,  a  50nm  thermally  deposited  film  was  also  scanned,  as  shown  in  Fig.  11b. 
The  evaporated  film  has  a  root  mean  square  (RMS)  roughness  of  3.27  nm  while  the 
epitaxially  grown  film  has  a  roughness  of  only  0.36  nm,  nearly  an  order  of  magnitude 
smalle.  Low-energy  electron  diffraction  (LEED)  and  reflection  high-energy  electron 
diffraction  (RFIEED)  patterns  shown  in  Fig.  11c  and  lid,  respectively,  taken  in-situ, 
further  confirm  the  long-range  single  crystalline  nature  of  these  films. 

We  performed  SE  measurements  and  analyses^^'^^,  to  multiple  epitaxial  and  thermal 
silver  films.  We  present  here  results  from  three  representative  films:  an  uncapped  40  nm 
epitaxially  grown  film,  a  45  nm  epitaxial  film  capped  with  1.5  nm  of  MgO  and  2  nm  of 
AI2O3,  and  an  uncapped  50  nm  thermal  film  deposited  at  a  rate  of  0.35  nm/s  serving  as 
control.  The  MgO/AkOs  cap  is  a  crucial  element  to  prevent  the  rapid  degradation  of 
epitaxial  silver  due  to  surface  oxidation  in  ambient  conditions,  and  we  confirmed  via  SE 
measurements  that  the  Ag  film’s  pristine  quality  persists  even  after  capping.  In  our 
fitting,  we  used  simple  Ag/Si  or  capping/Ag/Si  structural  models  for  the  epitaxial  films  as 
shown  in  Fig.  lie.  The  fitted  result  for £2  from  each  film  is  plotted  against  those  for  a 
40  nm  thermally  deposited  film  measured  by  JC  and  the  data  compiled  by  Palik  from 
multiple  samples.  In  this  work,  we  concentrate  on  the  energy  range  below  the  interband 
transition  threshold,  which  is  relevant  to  most  plasmonic  applications. 


In  the  energy  region  below  3.8  eV,  the  contribution  to  £2  mainly  comes  from  intraband 
transitions  typically  accompanied  by  scattering  between  electrons,  lattice  vibrations 
(phonons),  surface  roughness  and  grain  boundaries  inside  the  bulk.  Since  epitaxially 
grown  films  are  single  crystalline,  scattering  from  imperfections  is  expected  to  be  lower 
than  that  in  thermal  films.  Our  measurements  indeed  show  significantly  lower  loss  than 
JC’s  values  in  the  1.8-2.5eV  range,  while  the  control  film  shows  values  in  between  JC 
and  Palik  data.  The  residues  in  this  region  are  small  and  centered  around  zero  for  all 
films,  suggesting  that  our  model  fits  the  data  very  well.  In  the  lowest  energy  range  below 
1.5  eV,  the  error  becomes  large  due  to  reduced  detector  efficiency  and  our  extracted  £2 
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appears  to  be  larger  than  that  reported  by  JC.  We  note  that  the  original  JC  data  contain 
large  errors  in  this  same  energy  region  (gray  shadow  in  the  Fig.  11f-1h).  Thus  our 
measurements  are  consistent  with  JC’s  in  this  lowest  energy  region. 


The  residues  shown  in  Fig.  11f-1h  are  not  completely  random:  on  close  inspection,  in 
the  energy  range  near  3.7eV,  indicated  by  the  orange  circle  in  all  panels,  one  can  see  a 
sizable  peak  in  the  residue  for  the  thermal  film  (Fig.  11  h)  that  is  absent  in  both  epitaxial 
films  (Fig.  Ilf  and  g).  We  suggest  that  this  peak  is  associated  with  the  presence  of  grain 
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Fig. 11  .  AFM  scans  of  (a)  a  45nm  epitaxial  (2  nm  AI203/1 .5  nm  MgO  capped)  and  (b)  a 
50nm  thermally  deposited  (0.03  nm/s  deposition  rate)  Ag  film,  (c)  LEED  and  (d)  RHEED 
patterns  of  an  epitaxially  grown  Ag  film,  (e)  Layered  structures  of  our  Ag  film  samples  with 
and  without  oxide  capping.  Energy  dependence  of  £2  extracted  from  SE  measurements  on 
(f)  an  uncapped  40  nm  epitaxial  film,  (g)  a  45  nm  epitaxial  film  capped  with  1 .5  nm  of  MgO 
and  2  nm  of  AI203,  and  (h)  an  uncapped  50  nm  thermal  film  deposited  at  a  rate  of 
0.35  nm/s. 
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boundaries.  The  same  feature  around  the  same  energy  level  is  present  in  the  JC  data. 
Previous  theoretical  studies  modeling  the  effect  of  grain  boundaries  have  suggested 
that  they  lead  to  higher  loss  in  a  certain  wavelength  range  determined  by  the  average 
size^®.  These  earlier  calculations  have,  in  fact,  predicted  higher  loss  in  this  energy 
range. 
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Fig.  12.  Schematic,  experimental  results,  and  simulations  of  propagation  distance 
measurements,  (a)  Schematic  of  the  propagation  distance  setup,  (b)  SEM  image  of 
the  launching  and  output  slots  on  the  45  nm  thick  epitaxial  Ag  film,  (c)  Propagation 
measurements  for  two  excitation  wavelengths  (632  nm  and  880  nm).  Exponential 
curves  are  fitted  to  the  data.  The  inset,  plotted  on  semilog  scale,  compares  the 
simulated  propagation  distances  at  two  different  film  thicknesses  using  the  measured 
optical  constants.  The  green  dotted  line  marks  the  wavelength  at  632  nm.  (d) 
Simulation  of  the  mode  profile  and  the  SPP  excitation  at  the  interface  between  Ag  and 
the  oxide  capping  layers. 

Our  measured  optical  constants  suggest  improved  theoretical  limits  to  the 
performance  of  plasmonic  devices.  In  order  to  experimentally  demonstrate  these 
improvements,  we  measured  SPP  propagation  distances  over  the  45  nm  epitaxial  Ag 
film.  We  excited  and  detected  the  SPPs  in  reflection  geometry,  as  illustrated  in  Fig.  12a. 
Light  incident  from  an  oblique  angle  on  a  single  groove  launches  the  SPPs,  which  are 
subsequently  detected  at  a  series  of  output  coupling  slits  with  increasing  distance  from 
the  launching  site,  as  shown  in  the  scanning  electron  microscope  (SEM)  image  in  Fig. 
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12b.  We  used  two  different  incident  wavelengths  (632  nm  and  880  nm).  The  integrated 
optical  signals  from  the  output  grooves  are  plotted  as  a  function  of  propagation  length  in 
Fig.  12c.  The  experimental  data  were  fitted  with  simple  exponential  functions,  and  we 
extracted  intensity  propagation  distances  of  22  ±  5  fim  and  42  ±  3  for  SPP  at  632 
nm  and  880  nm,  respectively.  Analytical  calculations  predict  propagation  distances  of 
21  fj.m  at  632  nm  and  78  fim  at  880  nm  for  an  ideal  layer  with  the  newly  extracted 
permittivity,  as  shown  in  the  inset  of  Fig.  12c.  The  propagation  distance  measurement  at 
632  nm  is  in  excellent  agreement  with  the  calculated  value.  It  is  well  known  that  the 
SPP  propagation  distance  strongly  depends  on  the  film  thickness.  With  a  film  thickness 
of  200  nm,  for  instance,  the  predicted  intensity  propagation  distance  would  increase  to 
123  pm  at  632  nm  and  378  pm  at  880  nm,  respectively.  The  SPP  launching  mechanism 
and  expected  modal  profile  is  further  verified  by  conducting  full-wave  simulations  of  the 
geometry  (Fig.  12d),  employing  the  experimentally  retrieved  optical  parameters. 

We  suggest  that  future  theoretical  calculations  on  metamaterials  and  plasmonic 
devices  based  on  Ag  should  incorporate  the  new  optical  constants  reported  here,  as 
they  better  capture  the  intrinsic  properties  of  bulk  Ag.  Because  we  were  able  to  fit  the 
experimental  data  with  a  simple,  three-component  analytical  model,  these  newly 
extracted  optical  constants  will  facilitate  the  calculation  of  other  important  parameters 
such  as  the  material  Q-factor  and  group  velocity  in  Ag.  While  the  reported  optical 
constants  would  not  apply  to  low-quality  films  produced  using  thermal  evaporations,  Ag 
nanoparticles^^,  nanoshell^®,  and  nanoplatelets^^  synthesized  using  wet  chemical 
procedures  are  considered  to  be  of  single  crystalline  structure,  for  which  these  new 
optical  constant  values  are  expected  to  apply.  We  anticipate  that  these  high-quality 
epitaxial  silver  films  and  their  improved  optical  properties  will  have  a  significant  positive 
impact  on  the  fields  of  plasmonics  and  metamaterials  as  already  demonstrated  in  the 
case  of  nanolasers®°. 

C.6.  Progress  on  the  spin  wave  project 

We  have  mainly  three  accomplishments  in  the  spin  wave  project.  First,  we 
demonstrated  that  the  spin  wave  amplitude  can  be  controlled  by  a  direct  current  in  a 
bilayer  structure  composed  of  a  heavy  metal  and  a  magnetic  material  via  the  spin  hall 
effect.  Second,  we  demonstrated  the  use  of  the  micro-Brillouin  light  scattering  (^u-BLS) 
technique  as  a  local  temperature  sensor  for  magnons  in  a  permalloy  (Py)  thin  film  and 
phonons  in  the  glass  substrate.  Third,  we  demonstrate  that  the  spatial  decay  of  spin 
waves  launched  by  a  microwave  antenna  deviates  from  an  exponential  function  typically 
assumed  because  of  interference  effects.  We  focus  on  the  first  two  accomplishments  in 
this  report. 

C.6.1.  Spin  wave  amplification 
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Spin  waves  are  fundamental  excitations  in  magnetic  materials.  Understanding  the 
excitation,  propagation,  and  damping  of  spin  waves  is  critical  for  a  wide  range  of 
practical  applications  including  magnetic  memories,  microwave  oscillators,  and  logic 
devices  beyond  the  conventional  CMOS  technology.  For  example,  as  the  data  rate  in 
magnetic  memory  devices  reaches  gigahertz  range,  the  processional  dynamics 
associated  with  spin  waves  has  become  relevant  in  limiting  operational  speed  in  these 
devices.  SWs  may  also  serve  as  a  link  between  phase  locked  nano-oscillators  or 
carry  spin  current  in  magnetic  microstructures  subject  to  a  thermal  gradient  SWs 
have  also  been  proposed  for  use  in  an  information  bus  in  a  hybrid  logic  device  and  for 
building  logic  gates 

A  major  obstacle  to  many  of  these  applications  is  the  fast  damping  or  short 
propagation  length  of  SWs  in  metallic  magnetic  layers.  Current  induced  magnetization 
manipulation  via  spin  transfer  torque  (STT)  has  been  proposed  to  compensate  for  the 
damping  of  SWs  The  search  for  efficient  STT  materials  has  led  to  the  investigation 
of  metallic  bilayer  structures  consisting  of  a  magnetic  layer  and  a  nonmagnetic  layer 
with  strong  spin-orbit  coupling  (e.g.,  Pt  Ta  and  W  ).  In  the  nonmagnetic  layer, 
spin  polarized  current  is  generated  via  the  spin  Hall  effect  (SHE),  in  which  electrons  with 
different  spin  states  are  deflected  to  opposite  directions  Such  spin  polarized 

current  then  exerts  a  STT  on  the  adjacent  magnetic  layer,  enabling  control  of 
magnetization  dynamics  via  electric  control. 

We  investigated  the  electric  control  of  propagating  spin  waves  (SWs)  in  a  CoBFe/Ta 
bilayer  waveguide  as  illustrated  in  Fig.  13.  We  use  ^u-BLS  to  observe  the  amplitude  of 
the  SWs  in  the  bilayer  waveguide  and  its  change  due  to  an  applied  magnetic  field  and 
DCs.  After  removing  the  effect  of  heating  in  the  analysis,  we  observed  an  8%  change  of 
SW  amplitude  at  a  reasonable  DC  density  in  a  device  that  has  not  been  fully  optimized. 
Further  improvement  in  device  design  and  fabrication  may  eventually  lead  to  complete 
SW  damping  compensation  thus  opening  many  exciting  opportunities  in  spintronics 
and  magnonics.  This  accomplishment  resulted  from  collaborations  with  Prof.  Ralph  and 
Prof.  Burhman’  groups.  We  reported  this  result  in  Physical  Review  B  Rapid 
communication  in  2014. 
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Fig.  13:  Schematics  of  spin  polarized  current  generation  in  a  bilayer 
and  the  experimental  setup  of  micro-BLS. 

We  briefly  describe  our  experimental  approach.  The  iJ-BLS  technique  offers  high 
spatial  resolution,  sensitivity,  and  dynamical  range  as  demonstrated  in  previous 
experiments  The  BLS  signal  arises  from  the  inelastic  scattering  of  light  by  SWs, 
and  is  proportional  to  the  intensity  of  the  SWs,  i.e.,  proportional  to  the  SW  amplitude 
squared.  A  linearly  polarized,  single  frequency  laser  beam  at  532  nm  was  focused  to  a 
spot  size  of  ~  1  ijm  in  diameter  on  the  sample.  The  laser  power  was  approximately  1 
mW.  The  scattered  light  with  orthogonal  polarization  was  collected  and  sent  to  a 
Sandercock  multi-pass  tandem  interferometer  to  resolve  the  inelastic  scattering  from 
SWs. 

For  a  fixed  excitation  frequency  of  8  GHz  and  the  chosen  spatial  location  (shown  in 
Fig.  13),  we  scanned  the  applied  magnetic  field  and  DC  passing  through  the  bilayer 
waveguide.  We  plot  the  SW  amplitude  as  a  function  of  the  applied  magnetic  field  at  10 
and  -10  mA  in  Fig.  14(a).  There  is  a  clear  difference  in  the  SW  amplitudes  when  DCs  of 
opposite  sign  are  applied,  which  is  the  key  evidence  for  DC-controlled  SW  amplification 
and  attenuation.  The  direction  of  DC  for  the  SW  amplification  is  consistent  with  that 
expected  from  the  SHE,  as  illustrated  in  Fig.  13.  The  value  Hr  is  approximately  the 
magnetic  field  that  corresponds  to  the  peak  of  the  BLS  signal  in  magnetic  field 
dependent  spectra  (e.g.  Fig.  14(a)).  As  the  DC  increases,  there  is  a  clear  shift  in  Hr  as 
shown  in  Fig.  14(b),  consistent  with  the  heating  effect. 
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Fig.  14:  (a)  Magnetic  field  dependent  SW  amplitude  when  DCs  of  opposite 
directions  were  applied,  (b)  as  a  function  of  DCs.  (c)  Change  in  the 
propagating  SW  amplitude  at  a  particular  magnetic  field  as  a  function  of  DCs 
obtained  by  subtracting  SW  amplitudes  for  DCs  of  opposite  direction. 


To  isolate  the  change  in  the  SW  amplitude  due  to  the  STT  arising  from  the  SHE,  we 
remove  the  effect  of  heating  by  comparing  the  SW  amplitudes  at  two  DCs  of  the  same 
magnitude  but  opposite  directions  at  a  particular  magnetic  field  because  the  Joule 
heating  should  be  same  for  equal  magnitude  of  DCs  regardless  of  their  directions.  In 
this  context  we  define  >  where  A_(A+)  corresponds  to  the  SW  amplitude  for 

j  <0  {j  >  0).  The  percentage  change  in  SW  amplitude  from  AAat  H=  470  Oe  is  plotted 
in  Fig.  14(c).  An  approximate  linear  dependence  was  observed  up  to  the  maximal 
current  density  applied. 

We  now  estimate  the  change  in  damping  following  the  approach  outlined  in  The 
spin  wave  propagating  with  amplitude  it(x, t)  can  be  modeled  as  i^=|roo  + 

Vg  (^-i  ^  -  /cojj  u  -  iriu  -  iAr]u,  where  coq,  k^,  and  v  are  the  angular  frequency,  wave 

vector,  and  group  velocity,  respectively.  At]  represents  the  STT  induced  change  in  the 
SW  decay  rate  r?.  In  our  notation,  one  has  At]  <0  (At]  >  0)  for  j  <  0  (j  >  0).  The  SW 
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amplitude  at  x  =  L  when  a  DC  is  applied  can  be  written  as  Aj  -  Aoe  ,  where  i4o  is 
the  SW  amplitude  at  0  mA.  When  we  compare  the  cases  with  equal  DCs  of  opposite 

directions,  we  obtain  AA  =  tanh^^^. 

'’a 

From  the  data  presented  in  Fig.  14(c),  we  extracted  AA  at  different  DCs  and 
calculated  Aij  as  shown  in  Fig.  15.  In  this  calculation,  we  have  used  L  =  2  fim 
corresponding  to  the  distance  from  the  edge  of  the  antenna  to  the  detection  point.  The 
group  velocity  corresponding  to  the  SW  excited  at  8  GHz  was  calculated  from  the 
dispersion  curve.  One  can  also  calculate  the  change  in  damping  from  Aij  using  AaY{H  + 
ZnMgff)  -  Arj  (right  axis  of  Fig.  15).  If  we  use  an  approximate  value  of  a  =  0.008  from 
the  literature  the  maximal  change  in  damping  is  about  13%  at  /  =  10  mA, 

corresponding  to  the  current  density  of  6.3  x  10^  A/cm^. 
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Fig.  15:  The  change  in  decay  rate  (left  axis)  and  damping  constant 
(right  axis)  as  a  function  of  DC. 
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C.5.2  BLS  as  a  temperature  sensor 

We  demonstrate  that  the  BLS  spectra  can  be  used  as  local  temperature  sensors 
for  magnons  and  phonons.  Such  local  temperature  sensors  may  be  applied  to 
investigate  spin  caloritronic  and  thermal  transport  phenomena.  For  example,  a  key 
assumption  in  the  current  understanding  of  thermally  driven  spin  currents  is  that 
magnons  temperature  may  not  reach  thermal  equilibrium  with  the  phonons  locally.  This 
difference  between  the  phonon  and  magnon  temperatures  underlies  the  so  called 
“phonon-magnon  drag”  that  leads  to  the  directional  magnon  transport  or  spin  current  in 
magnetic  insulators  and  semiconductors^°^'^°^.  Direct  experimental  confirmation  of  this 
phenomenon  is  currently  lacking  most  likely  due  to  the  small  temperature  difference 
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between  these  two  temperatures  in  a  sample  with  a  modest  temperature  gradient^°^ . 
This  accomplishment  resulted  from  the  close  collaboration  between  Prof.  Li  Shi  and  is 
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Fig.  16.  Experimental  schematics,  (a)  Uniform  heating,  (b)  A  thermal  gradient  is  applied 
laterally  to  the  substrate  and  perpendicular  to  the  magnetic  field. 

published  in  Applied  Physics  Letters. 

Our  experimental  setup  is  illustrated  in  Fig.  16.  The  sample  consists  of  a  Py  film 
sputtered  on  a  glass  substrate.  The  film  is  nominally  60  nm  thick  and  coated  with  10  nm 
of  SiOx  to  prevent  oxidation.  The  substrate  is  a  rectangular  piece  of  glass  approximately 
1  mm  wide  and  3  mm  long.  A  magnetic  field  of  500  Oe  was  applied  parallel  to  the  short 
side  of  the  sample.  For  the  uniform  heating  measurements  (Fig.  17(a)),  the  sample  was 
mounted  on  a  Peltier  heater  that  was  much  larger  than  the  sample.  For  the  thermal 
gradient  measurements  (Fig.  17(b)),  one  end  of  the  film  was  mounted  on  a  resistive 
heater  while  the  other  end  was  mounted  on  a  heat  sink.  A  thermal  gradient  was  applied 
along  the  length  of  the  sample  and  perpendicular  to  the  external  magnetic  field. 

We  first  calibrated  how  the  magnon  frequencies  shift  as  a  function  of  temperature 
when  the  sample  is  heated  uniformly,  as  shown  in  Fig.  17.  The  spectra  are  dominated 
by  two  strong  peaks  corresponding  to  the  first  two  perpendicular  standing  spin  wave 
(PSSW)  modes  (n  =  1,  n  =  2)  of  the  Py  film.  The  PSSW  modes  are  characterized  by 
sinusoidal  amplitude  profiles  in  the  film  thickness  direction  with  the  number  of  nodes 
corresponding  to  n.  For  the  relatively  small  in-plane  wave  vectors  probed  in  our 
experiments  with  normal  incident  light,  the  frequencies  of  the  PSSW  modes  can  be 
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calculated  according  toJ°  ^  (y)  ) (^ext  +  ^  (y)  +  ,  where  y  is 

the  gyromagnetic  ratio  (2.93  MHz/Oe),  A  (1.0  erg/cm  at  -300  is  the  exchange 
constant,  Next  is  the  applied  external  field,  L  is  the  thickness  of  the  ferromagnetic  film, 
and  Ms  is  the  saturation  magnetization.  The  temperature  dependence  of  the  exchange 
constant  and  saturation  magnetization  cause  the  mode  frequency  to  vary  with 
temperature. 

The  magnon  frequency  was  extracted  by  fitting  each  PSSW  peak  in  the  BLS 
intensity  with  a  Lorentzian  lineshape,  as  for  the  magnetic  field  dependent 
measurements.  The  measurements  were  fit  individually  and  the  centers  of  the  peaks 
were  averaged  for  24  independent  measurements.  The  error  in  frequency  quoted  for 
each  temperature  is  the  standard  deviation  of  the  mean.  Over  the  temperature  range 
studied,  the  frequency  decreases  linearly  with  increasing  temperature  due  to  the 
reduced  saturation  magnetization  and  exchange  constant  at  higher  temperatures. 
Linear  fits  to  the  data  yield  slopes  of  (-8.4  ±0.4)x10'^  GHz/K  and  (-16.1  ±0.8)x10'^ 
GHz/K  for  the  n  =  1  and  n  =  2  modes,  respectively.  The  average  error  in  the  frequencies 
for  the  two  modes  are  0.02  GHz  for  n  =  1  and  0.05  GHz  for  n  =  2.  Therefore,  the 
measured  temperature  sensitivities  are  about  2.4  K  and  3.1  K  for  the  n  =  ^  and  n  =  2 
modes,  respectively.  The  temperature  sensitivities  quoted  in  this  study  can  be  further 
improved  simply  by  extending  the  measurement  time,  but  the  measured  values  already 
compare  favorably  with  other  nanoscale  thermometry  techniques^°^. 
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Fig.  17:  Temperature  dependent  thermal  magnon  spectra  in  a  uniformly  heated  film,  (a) 
Thermal  magnon  spectra  with  increasing  temperature  from  the  bottom  to  the  top.  Extracted 
magnon  frequencies  as  a  function  of  temperature  for  the  (b)  n  =  1  and  (c)  n  =  2  modes.  The 
solid  lines  are  linear  fits,  which  allow  us  to  determine  the  frequency  shift  as  a  function  of 
temperature. 


We  now  investigate  the  feasibility  of  using  the  BLS  spectra  as  a  local  magnon 
temperature  sensor.  To  do  so,  we  establish  an  in-plane  temperature  gradient  by 
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mounting  the  sample  to  a  heater  and  heat  sink,  as  illustrated  in  Fig.  17(b).  We  measure 
the  magnon  spectra  at  different  locations  along  the  temperature  gradient  by  moving  the 
sample  on  a  motorized  stage  in  10  steps  of  160  pm  each.  We  performed  24 
independent  measurements  at  each  position,  averaged  the  peak  positions  fitted  with  a 
Lorentzian  lineshape,  and  quoted  the  standard  deviation  of  the  peak  positions  as  the 
error  bar.  To  account  for  the  spatial  variation  of  the  magnon  frequencies  due  to  spatially 
non-uniform  properties  of  the  film,  we  subtracted  the  measured  magnon  frequencies  at 
each  location  without  the  thermal  gradient  from  those  with  the  thermal  gradient.  The 
resulting  magnon  frequency  shifts  for  the  n  =  2  and  n  =  ^  PSSW  modes  are  displayed  in 
Fig.  18(a)  and  18(b),  respectively.  The  spatial  dependence  of  the  magnon  spectra 
corresponds  to  a  linear  temperature  profile,  as  expected  in  this  simple  geometry.  We  fit 


Fig.  18:  The  extracted  magnon  frequency  shift  for  the  (a)  n  =  1  and  (b)  n  =  2  modes  as  the 
beam  was  scanned  along  the  temperature  gradient.  The  magnon  frequency  without  a 
thermal  gradient  is  subtracted  from  that  with  the  thermal  gradient.  The  approximate  linear 
shift  in  magnon  frequency  is  consistent  with  a  linearly  varying  temperature  profile 

the  frequency  shifts  to  a  linear  function  and  obtain  shifts  of  (-8.9  ±1.3)x10'^  GFIz/mm  for 
the  n  =  1  mode  and  (-20  ±2)xl0'^  GFIz/mm  for  the  n  =  2  mode.  Using  the  uniform 
heating  data,  the  measured  frequency  shifts  correspond  to  a  -10.6  ±1.2  K/mm 
temperature  gradient  between  the  hot  and  the  cold  ends  for  the  n  =  1  mode  and  -12.6 
±1 .2  K/mm  as  measured  with  the  n  =  2  mode.  We  note  that  the  temperature  gradients 
extracted  from  the  two  PSSW  modes  are  consistent. 

To  complement  the  magnon  temperature  measurement,  we  measure  the 
acoustic  phonon  spectra  of  the  glass  substrate.  The  spectrum  contains  one  peak 
corresponding  to  the  frequency  of  the  acoustic  phonons  in  glass  with  wave  vectors  of 
4-nN/K,  where  A  is  the  laser  wavelength  and  N  is  the  index  of  refraction  at  the  laser 
wavelength.  Assuming  N  =  1.5,  the  measured  frequency  of  35  GHz  implies  a  speed  of 
sound  of  6  km/s,  which  is  quite  reasonable  for  glass.  We  first  calibrated  the  change  in 
the  phonon  spectra  as  a  function  of  temperature  under  uniform  heating  conditions  with 
the  setup  shown  in  Fig.  17(a).  As  the  temperature  was  raised,  a  systematic  increase  in 
the  intensity  of  the  phonon  peak  was  observed  (Fig.  19(a)),  although  the  peak  position 
did  not  shift  in  frequency.  Because  the  ratio  between  the  phonon  energy  and  kT  is 
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Fig.  19:  (a)  Phonon  spectra  of  the  glass  substrate  collected  at  three  different 
temperatures  under  uniform  heating,  which  shows  how  the  peak  intensity  increases  as 
the  temperature  increases,  (b)  Percent  change  of  the  intensity  of  the  phonon  peak  as  the 
temperature  was  increased  under  uniform  heating,  (c)  Percent  change  of  the  intensity  as 
the  beam  was  scanned  along  the  temperature  gradient,  which  shows  a  linearly  varying 
temperature  profile. 

small,  the  Bose-Einstein  distribution  for  the  phonon  mode  is  reduced  to  the  classical 
limit  that  depends  linearly  on  temperature,  i.e.,  <n>  =1/(exp(/76u//fT)  -1)  «  kVfiuj. 
Consequently,  the  intensity  of  the  phonon  peak  in  the  BLS  spectra  depends  linearly  on 
temperature,  as  shown  in  Fig.  19  (b).  The  intensities  plotted  in  Fig.  19(b)  are  the 
averaged  value  from  24  measurements  with  the  standard  deviation  displayed  as  the 
error  bar.  A  linear  fit  yields  that  the  phonon  intensity  increased  by  22  ±2  %  as  the 
temperature  increased  from  296  K  to  360  K,  also  an  increase  of  22  %.  The  slope  of  the 
linear  fit  is  0.35  ±0.03  %/K.  To  measure  the  local  temperature  of  the  sample  with  a 
thermal  gradient,  we  subtracted  the  intensities  at  each  point  with  and  without  the 
thermal  gradient  to  compensate  for  non-uniformities  in  the  substrate.  The  data  shown  in 
Fig.  19  (c)  is  fit  with  a  linear  function  with  a  slope  of  6.3  ±1.1  %/mm.  Using  the  slope 
obtained  from  the  uniform  heating  data,  this  change  in  BLS  intensity  of  the  phonon 
mode  implies  a  temperature  gradient  of  18  ±3  K/mm  over  the  region  measured,  which  is 
somewhat  larger  than  the  value  extracted  from  the  magnon  spectra  but  still  within  two 
standard  deviations.  Further  work  is  required  to  understand  whether  a  difference  exists 
between  the  phonon  and  the  magnon  temperature  gradients  due  to  local  non¬ 
equilibrium  between  phonons  in  the  glass  substrate  and  magnons  in  the  Py  film. 

C.6.3  Deviation  from  exponential  decay  of  spin  wave 

We  investigated  the  propagation  of  magnetostatic  surface  spin  waves  (MSSW)  in 
a  continuous  Py  thin  film  excited  by  an  asymmetric  coplanar  waveguide  antenna.  The 
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spin  waves  were  excited  in  the  Damon-Eshbach^°®  configuration,  wherein  a  static 
magnetic  field  was  applied  parallel  to  the  antenna  and  perpendicular  to  the  propagation 
direction  of  the  spin  waves.  While  mapping  the  spatial  dependence  of  the  surface  spin 
wave  intensity,  we  observed  a  significant  deviation  from  the  simple  exponential  decay 
due  to  damping:  an  oscillatory  signal  is  superimposed  on  top  of  the  exponential  decay. 
To  model  the  observed  complicated  spatial  dependent  spin  wave  intensity,  we  use  a 
closed-form  expression  that  includes  contributions  from  the  propagating  spin  wave  and 
a  background  magnetization  with  spatially  uniform  phase.  This  simplified  closed-form 
model  allows  us  to  fit  the  experimental  data  and  to  extract  several  key  parameters 
including  the  wave  vectors,  the  group  velocities,  and  the  propagation  lengths  of  the 
excited  spin  waves.  This  result  is  published  in  Applied  Physics  Letters. 
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